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Abstract. We consider stability and approximate reconstruction of Rie- 
mannian manifold when the finite number of eigenvalues of the Laplace- 
Beltrami operator and the boundary values of the corresponding eigenfunc- 
tions are given. The reconstruction can be done in stable way when manifold 
is a priori known to satisfy natural geometrical conditions related to curva- 
ture and other invariant quantities. 

I. Introduction In this paper we consider the questions of stability and 
approximate reconstruction in the inverse boundary spectral problem which 



is also called the generalized Gelfand inverse problem [13| for Riemannian 
manifolds. To formulate the problem and the main results, we need to 
introduce some basic notations. We will denote by (M, g) an unknown, 
m-dimensional, compact connected Riemannian manifold with a (smooth) 
metric g and non-empty boundary dM. Our goal is to find {M,g) from 
boundary data. The boundary dM is itself an (m — l)-dimensional compact 
differ entiable manifold. We note that we do not assume the knowledge of 
the metric i*{g) generated on dM by the metric g , where i : dM M is 
an embedding or even the corresponding area element dSg. 

Because in the inverse problems the boundary dM is usually known, 
we will consider a class M = M^jv/ of compact, connected Riemannian 
manifolds which have the same, i.e. diffeomorphic, boundary, dM. 

Let — be the Laplace operator on (M, g) with Neumann boundary 
condition. Denote by (fj; j = 1,2,...} the complete set of eigen- 
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values and corresponding orthonormal eigenfunctions of — A^, = Ai < 
A2, ...,ifi= vol~^/^(M, (7), where vol stands for the volume of {M,g). 

Denote by B the set of sequences {nj, tpj] j = 1, 2, . . .} where fij G R 
and ipj € L^{dM) and by P : M ^ B, the map 

V{M,g) = {Xj, ifjlQM, j = 1,2,...}. 



Definition 1 The collection {Xj, ipj\dM', j = 1) 2, . . .} is called the boundary 
spectral data of (M, g) and the map T> - the boundary spectral map. 

It was shown by Belishev-Kurylev 0] who used the boundary control 
method Q and the unique continuation result of Tataru |35| that 



Theorem 1 The map 2? : M ^ B is injective. 

(For earlier uniqueness results in the Gelfand inverse problem for isotropic 
operators obtained by the boundary control method see, e.g. p). The com- 
plex geometric optics method was applied to this problem in |2^], p^].) 

In this paper, we will first analyse the question of stability of the inverse 
problem, i.e. the question of continuity of'D~^. Later, we will also describe 
a procedure of an approximate reconstruction of {M,g). 

Alessandrini Q , [|| considered an inverse boundary problem for a Schrod- 
inger operator, —A + q in a bounded domain of R™' and obtained a log-type 
stability estimate in the Gel'fand inverse problem (see also p3| ). 

Alessandrini and Sylvester |^] and Stefanov and Uhlmann considered 
inverse boundary value problems for the wave equation which are closely 
related to the Gel'fand inverse problem. Let 

utt + a{x,D)u = 0, in 17x[0,T], u\t<o<0, ^i|anx[o,T] = /, 

where a{x, D) = —A + q in and a(x, D) = —dia^Wj in [^] and the inverse 
data is given as a non-stationary Dirichlet-Neumann map. For sufficiently 



large T and, in |3l[| a*-' close to 6'^^ , they showed Holder-type stability. These 
results raise the question of an optimal stability estimates in the Gel'fand 
inverse problem and inverse boundary value problem for the wave equation. 
Although exact stability estimates lie outside the scope of this paper which 
is devoted to the analysis of geometric conditions upon (M, g) to provide 
continuity of we believe that Holder-type estimates could be valid only 
under rather strong requirements upon (M, g), like strong geodesic property 



(see e.g. [|^], |^^). Further results on stability for different types of inverse 
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boundary value problems could be found in [16|. For counterexamples see 



e.g. HI 



Clearly, to discuss the continuity of it is necessary to define appro- 
priate topologies in M and B. However, it is well-known (see e.g. ||Tl|| ) 
that inverse problems are, in general, ill-posed. Therefore, to gain stability 
we need to know a priori that the unknown object, which we are going to 
reconstruct, belongs to some compact class K (see e.g. To obtain 

stability results, the topologies of M and B have to satisfy two properties. 
First, the map 

V :K^B 

has to be continuous. Second, D has to be injective. When these two 
properties are satisfied and K C M is a compact set then the continuity of 
the map 

: V(K) K 

follows from basic results of general topology. In the first part of the paper 
we carry out the constructions of the appropriate topologies. We would like 
to note that a proper topology on M is the Gromov-Hausdorff topology 
and the (pre)-compact subsets we actually use are the manifolds of bounded 
geometry (see e.g. |14], [^], [29|). The topology on B is defined from 



the boundary spectral convergence. 

In the second part of the paper our attention will be focused on the 
problem of an approximate reconstruction of (M, g) when we do not know 
all boundary spectral data. Rather, we know only the first eigenvalues, Xj < 
5^^ with some small 6 > and boundary values of the corresponding 
eigenfunctions and, furthermore, we know them with some error. Namely, 
we have a finite collection of pairs 

{fij,iPj; j = 1, . . . , n{6-^)}, fij G R, Vi G L\dM) 

and know that they are close to {Xj, ^PjldM', j = ■ ■ ■ , f^i^^^)}, where 
n{5'^) is the number of eigenvalues that are smaller than 5^^. Our defini- 
tion of topology on B is adjusted to make the intuitive notion of closedness 
of to {Xj, ifjldM] 3 < rigorous. 

In the second part of the paper we describe the stable reconstruction 
procedure but omit the proof showing that this procedure is stable. These 
proof are to be published elsewhere. We note that to carry out the recon- 
struction procedure we will need to work under stronger requirements on 
(M, g) than those that are needed for stability. Similarly, in the first part 
we will only formulate necessary results omitting all proofs. The proofs 
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and detailed description of the reconstruction procedure will be given in a 
forthcoming paper. 

II. Main Results. Stability. We start with a proper topology on M. 

Definition 2 (Gromov-Hausdorjf topology). Let £ > 0. The Riemannian 
manifolds {M^,g'^) G M, i = 1,2 are e-close in the Gromov-Hausdorff topol- 
ogy, 

dGHm\g'), {M^g^))<e, 
if there are e-nets {x* ; j = 1, . . . , J(e)} C M*, i = 1, 2, such that 

\dg^ (xj, x^) — dg^ix'^j , x\)\ <£, j,k = 1, . . . , J{£), 

where dgi{-, ■) stands for the distance on (M*,^'). 

We remind that {xj; j = 1, . . . , J(e)} C M is an e-net if for every x G M 
there is Xj such that dg{xj,x) < e. 

Remark 1. In the future we will identify diffeomorphic manifolds assuming 
implicitly that desired statements are valid after an automorphism of M. 

Definition 3 (Riemannian manifolds of bounded geometry). For any A, 
D, if) > 0, M(A, D,zo) C M consists of Riemannian manifolds {M,g) G M 
such that 

i) \\VRm{M,g)\\ + \\Rm{M,g)\\ < A, Hi) diam{M,g) < D, 

ii) \\VS{M,g)\\ + \\S{M,g)\\ < A, iv) inj{M,g) > i^. 

Here Rm is the Riemannian curvature tensor on (M,g) considered as 
a 4-linear form on TM, V is the covariant derivative, and S is the second 
fundamental form with respect to the inner product, i.e. a quadratic form on 
TdM. Furthermore, diam is the diameter of (M, g) and inj is the injectivity 
radius of {M,g). Condition iv) consists of three different conditions: 

(a) Let X G dM. Denote, as usual, by BgM{x,r) the open metric ball (on 
dM) of the radius r with center in x. Then, for any r < io, BQM{x,r) is a 
domain of normal coordinates on dM centered at x. 

(P) Let X G M, d{x,dM) > io. Then, for any r < io , BM{x,r) is a domain 
of normal coordinates on M centered at x. 

(7) Let X G M, d(x, dM) < io. Then, for any r < io, the cylinder C(x, r) is a 
domain of boundary normal coordinates. Here C{x, r) consists of all points 
y £ M such that d{y, dM) < r and the unique boundary point z{y) G dM, 
nearest to y is in the ball BQMiz{x),r). 

Next we introduce topology on the set B. 
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Definition 4 (Boundary spectral topology.) Let 6 > 0. Collections 

Vj; J = 1,2,...} GB, i = l,2 
are 5-close if there are disjoint intervals 

Ip = {ap, bp) C (-(5-1 -5, 5-^ + 6), p=l,..., P, 

such that 

i) bp — ap < 5, 

ii) For any fij, i = 1,2 with /i*- < 6^^ there is p such that ^* G Ip. 
Hi) For any p, the total number of fi^j inside Ip coincide, i.e. 

= nl{= Up), 

iv) For any p there is a unitary matrix Up = [u^^^^-^ such that 

Here, ^'^ is the vector-function ' ' "^^^n ))' ^^^^ < ••• < j{np), with 

P'j(k) £ Ip for k = 1, . . .,np. 

Clearly, condition iv) depends on the choice of the boundary measure 
used in the definition of L^{dM). However, due to the compactness of dM 
different smooth measures on dM determine equivalent L^-norms. 

We note that such topology was introduced by Alessandrini [|^, [^] who 
studied stability in the Gel'fand inverse problem for a Schrodinger operator. 

We are now able to formulate the principal stability result for inverse 
boundary spectral problem. 

Theorem 2 For any A, D, iq > 0, the map exists and is continuous on 
M{A,D,io). That is, there is 6 > such that if{M\g^) £ M(A, io), i = 
1,2, and their boundary spectral data, 

{a], 9?]|sm; j = 1,2, . . .} and {A|, (^||aAf; j = 1, 2, . . .} , 

are 5-close then and are diffeomorphic. Moreover, for any a G (0, 1), 
— > g^ in C^'" -topology when 5^0. 
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Remark 2. The convergence of the metrics in the C^'"-topology means 
that in a proper coordinate system on M = Mi, i = 1,2, the metric tensors 
gij converge in the C"^'"-sense, see e.g. [p^]. 

Remark 3. It follows from definition |^ that to apply Theorem Q it is 
sufficient to know only the parts of the boundary spectral data of (M^ , ) 
and (M^,g'^) that correspond to the eigenvalues A*- < 6~^. 

Theorem Q and the fact that in the class M(A, D, ig) the Gromov- 
Hausdorff topology is equivalent to the Lipschitz topology (see ([^) implies 
the following corollary that describes the Lipschitz convergence of Rieman- 
nian manifolds. 

Corollary 1 Let {M\g'-) £ M(A, L>, zq), i = 1,2. Then, for any e > 
there is 6 > such that if the boundary spectral data of these manifolds are 
5-close then and are diffeomorphic and 

, 1 d„i(x,y) 
l + g -i< 9 < 1 + 6 , x,yGM. 

dg2{x,y) 

Remark 4. It should be noted that M(A, Z),io) is not closed in M in the 
Gromov-Hausdorff topology. Therefore, to obtain the desired stability re- 
sult we should first extend Theorem |l] onto the closure M(A,Z),zo)- The 
set M(A, D, io) consists of some class of C^'"-smooth Riemannian manifolds. 
Although the boundary control method that was instrumental to prove The- 
orem ^ is, in general, not applicable to C^'"^ Riemannian manifolds it was 
possible to show that 

V : M{A,D,io) B 

is injective (see section IV). 

In Theorem ^, the condition iii) in M(A, D, io) is automatically satisfied 
from an estimate of A2 in [^]. The other conditions are natural for the 
stability of the inverse problem because, otherwise, the convergence of the 
boundary spectral data does not imply the convergence of the corresponding 
Riemannian manifolds. We illustrate this thesis with the following examples. 

Example 1. Let {M,g) be a two-dimensional manifold obtained by the 
following surgeries of the two-sphere S'^. Add a handle H near the North 
pole and cut-off a small piece B to make a boundary near the South pole. 
If the metric size of the handle tend to uniformly then, in the bound- 
ary spectral topology, the boundary spectral data of (M, g) tend to the 
boundary spectral data of the sphere with a hole but without a handle (see 



e.g. [34|). This shows the necessities of i) and iv) in 'M.{A, D,i()) 
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Figure 1: Two representations of the same manifold: Cylinder in R'^ pinched 
in the middle and its representation as an annulus in with variable metric. 
When e ^ the arclength of circles in the region A tend to so that the 
manifold splits into two components. 

Example 2. Let {M,g) be cylinder in R'^ with the boundary made of 
two components, dM = dMi n dM2 (see Fig.l). We pinch the cylinder in 
the middle so that the diameter tends to 0. Therefore, the manifold even- 
tually "splits" into two disjoint semispheres. Meanwhile, in the boundary 
spectral topology on dM, the boundary spectral data of (M, tend to the 
boundary spectral data of this disconnected manifold. Again, conclusion of 
the theorem is no more valid. 

This second example may be also described in terms of a perturbation 
of the metric in an annulus A{1, 4) = {(r, ^) : 1 < r < 4} , where (r, 6)) are 
the polar coordinates in R^. We define the length element by 

2 

die = dr^ + ■■ ^ , , d9^ x{r) e C^, 

1 + e '^xin 

where 

X = for r < 2 or r > 3, x > for r G (2, 3). 

Then, for £ ^ each circle r G (2, 3) shrinks to a point so that the internal 
part near r = 1 of A becomes separated from the external part near r = 4. 

We would note that although this example does not precisely satisfies 
the definition of 'M.qm because the limiting manifold is no longer connected. 
However, this definition can be extended onto the case of, probably discon- 
nected Riemannian manifolds with the same boundary dM which have no 
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components without the boundary if we change the condition of bounded- 
ness of diam(M) to the boundedness of vol{M). 

Example 3 Take X and Y to be isospectral manifolds without boundary of 
different topological types (see e.g. jS^) and Z he a Riemannian manifold 
with boundary dZ. Let be a Riemannian manifold obtained by connect- 
ing Z with X by a thin tube of radius of length e. Similarly, let Ngr be 
a Riemannian manifold obtained by connecting Z with y by a thin tube of 
radius of length e. When e — > 0, the boundary spectral data of Me tend 
to the union of the boundary spectral data of Z and {/Xj,0}, where is 
the spectrum of X. Clearly, the boundary spectral data of tend to the 
same limit. Thus, for any 6 > 0, there exists > such that the boundary 
spectral data of and Ne are (5-close when e < eo- However, and 
are of different topological types. 

In these examples ||-R?n'|| oo and inj when e ^ which show the 
importance of conditions i) and iv) of the theorem. Similar examples can 
be given to illustrate the importance of condition ii). 

III. Main results. Approximate reconstruction. Definition ^ of 
the Gromov-Hausdorff topology in M can be easily extended to the set of all 
compact metric spaces {X, d). Thus, two Riemannian manifolds are 
e-close in the Gromov-Hausdorff metric if they possess e-nets X* which, be- 
ing considered as finite metric spaces with distance inherited from (M^,g'^), 
are e-close in the Gromov-Hausdorff topology. Therefore, to approximately 
reconstruct the Riemannian manifold (M, g) we should construct a metric 
space (y, d),Y = {yi, . . . , yj} which is e-close to (M, g), i.e. 

dGHm,9), (y,d))<e. 

To construct a metric space approximation to {M,g) we will assume 
further regularity properties of the class of manifolds. 

Theorem 3 Let A, D, iq > 0. Then for any e > there is 5 > such that 
if {/ij, tpjldM'^ J = 1)2, ...} G B is 5-close to the boundary spectral data 
{Xj, (pjldAf, j = 1,2, . . .} of a Riemannian manifold, [M,g) G M(A,D,io) ; 
then there is a finite metric space {Y, d) such that 

dGHm,g), {Y,d))<£ 

Moreover, there is a constructive algorithm to find (Y,d) from {/x*-, -0}; J = 
1,2,. ..}GB. 
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The algorithm to construct (Y, d) is described in section V. 

Remark 5. By definition ^ the 5-closedness does not involve {/i*-, -0]} with 
> 6^^. Similarly, the construction of {Y,d) does not use these data. 

IV. Geometric convergence and stability. The proof of Theorem ^ 
is based upon the following result. 

Proposition 1 For any A,D,iQ > 0, the set M(A, D,io) is pre-compact in 
the Gromov-Hausdorff topology. Its closure, M(A, D, iq) consists of differ- 
entiable manifolds with C^^'^-smooth metric, where a G (0,1) is arbitrary. 

If a sequence of Riemannian manifolds {M"',g'^) converges to {M,g) in 
the Gromov-Hausdorff topology on M(A, D, io), there is no > 1 such that 

i) For n > no, there is a diffeormorphism Fn : M ^ M", i.e. M" and M 
are diffeomorphic. 

ii) for any a E (0, 1), F*{g^) converge to g in the C^'"-topology. 

Remark 6. Proposition |^ is obtained by applying the ideas developed 
by M. Anderson (Q, see also |jl5|) in the interior of the manifold and the 
fundamental equation of Riemannian geometry (see [^]) near the bound- 
ary. The latter equations is also known as the Riccati equation for the 
second fundamental form along normal geodesies. This equation gives the 
desired regularity estimates for the metric tensor in the boundary normal 
coordinates. 

The geometric convergence described by Proposition [l| yields the conti- 
nuity of the direct problem. 



Proposition 2 T) : M(A, D,io) ^ B zs continuous. 



Remark 7. Proposition [T| was proven by Kodani ||2T| in the case of weaker 
regularity. Namely, Kodani showed that the manifolds satisfying 

i) \\Rm{M,g)\\ < A, iii) diam(M,c/) < D, 

ii) ||5(M,5)|| < A, iv) inj(M,ff) >io. 

are pre-compact in the Lipschitz topology rather than the C^'"-topology. 
This result is sufficient for continuity of T>. This can be obtained by means 
of the perturbation theory of quadratic form (see e.g. 1 12] in the case of man- 
ifolds without boundary). Therefore, the result of Kodani would have been 
sufficient to prove the continuity of 1)-^ : Vi MjA, D,!^) M(A,L',io) if 
it were possible to show the injectivity of V on M(A, D,io). Unfortunately, 
the method used to prove Theorem |l] fails in this case in its parts related to 
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the approximate controllability result by Tataru p5| | and also to the injec- 
tivity of the boundary distance map (see the next section for the definition 
of this map). 

Proposition |l], on the contrary, guarantees a stronger regularity of the 
Riemannian manifolds in M(A, D, iq). It was used to extend Theorem |^ on 
the class 

M= (J M{A,D,io). 

A,D,io 

Theorem 4 The map I? : M ^ B is injective. 

V. Construction of a finite e-net. The procedure of constructing an 
approximation (Y, d) consists of several steps which we will explain sepa- 
rately. We will also point out smoothness requirements which are necessary 
to carry out different steps of the procedure. 

a. Construction of Fourier coefficients of waves. For simplicity, we 
assume here that the metric on the boundary, i*{g), i : dM — > M is given. 
Consider an initial-boundary value problem associated with the Neumann 
Laplace operator, 

u{t - AgU^ = 0, in M X [0, D], 

uf\t=Q = u{\t=Q = 0, (1) 

where Cp([0,D], L^{dM)) consists of piecewise C^-functions oft. We intro- 
duce the wave operator given by 

W\f) = uf{;t)eCH[0,D],L\M)) 
and, if we denote by the A;-th Fourier coefficient of u-^ , 

u{{t) = {uf{t), = f j si{y,t')f{y,t')dt'dSg{y), (2) 

Jo JdM 

where dSg is the boundary area element and 

t( ,/N sin(\A^(t - t')) / N 



10 



This formula makes possible to find approximately the first Fourier coeffi- 
cients of the waves u^^ if we know {fij^t/jj; j = 1, . . . , n{6^^)}. 

b. Construction of domains of influence. Let F C dM be open. 

Consider the waves n-^ (r) with / G C^(F x [0, r]), r > 0. Clearly, 

supp[u-^ (r)] C M(F, t) ={x£ M : dg{x, F) < r}. 



By the fundamental result of Tataru see also Ch.2.5), we have 

Theorem 5 Assume that the metric tensor g is -smooth. For any r > 
and F C dM, 

cIlHm) {u^ir) : / e C^iF X [0,r])} = L2(M(F,r)), 

where C L'^{M) consists of functions with support in O. 

Let r] > 0, D > diam(M) as in part iii) of Definition ^ and F;, / = 1, . . . , L, 
be open subsets of dM satisfying 

diam(FO < 7?, F; n Ffc = 0, dM = uT7. 

Let 

a = {ai,- ■ ■ ,aL), ai<D/r], (3) 
be a multi-index. Denote by Sq, C dM x [0, D] the set 

= U (F; X - aiTj, D]) n (9M X [0, D]) 
I 

and by M^ C M the subdomain 

M^ = \jM{Fi,air]). (4) 

By Theorem ^, for any cr > there is / = /q G C(f^(Sa) such that 

where xn is the characteristic function of a set $7. 

It then follows from the general results of control theory that there is a 
finite linear combination /, 

J 

f{x,t) =Y,ajsf{x,t)xi:Ax,t), xedM, (5) 
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such that 



IWiD) - XMc,Vi\\l^m) < (6) 

The coefficients Uj of the function (^) which satisfies equation (^) can be 
found using the boundary spectral data {A^, fk\dM] k = 1, . . . , K}, K > J 
by means of some variational procedure. 

To define this procedure rigorously, let HN^^a) C C^{[0,D], L^{dM)) 
be an A^-dimensional subspace spanned by the functions {s^{y, t) xt,^ iu^ 
n=l,...,N}. 

Let Anj^kH) be a functional on H]\f(I!,a), 

Anj,kU) = (7) 
= snp{\\{PKWf-XM^^i,Wh)\\ : h^Hii^^), WPx^hW < l} . 
Here P/v is the orthoprojector in L?'{M), 

Pnu = ^ {u,ipn) Vn- 

n<N 

Clearly, due to (P) the functional j4(/) can be evaluated in terms of {Aj, ^Pj\dM] 
j = 1, . . . , if n{6^^) > max(A^,/,K). As we have {^j^ijjj]] = 

1, . . . , n((5~^)} rather then {Aj, ^j\dM'-, 3 = 1, ■ • • > ^(f^""^)}, the functional 
A{f) can be evaluated only approximately. 



Theorem 6 For any o" > there are parameters C, N,I,K,6, N < I < 
K < n[6^^), such that if {fij, ipj} is 5-close to the boundary spectral data 
{Aj, ^j\dM} of a Riemannian manifold {M,g), then a minimizer f* G H]\f 
of the functional Aj^j^k such that 

ll/*llL2(gMx[0,D]) < C, 

satisfies the equation 

\\Wf* -XM^Vi\\<a. 



By Proposition || the boundary spectral data depend continuously on 
the Gromov-Hausdorff distance restricted on M(A, D,io)- Consider 



w 



as a function of {M,g) E M(A, D, io). By Proposition ||, the perturbation 



theory |1£] shows that the wave u-^j^j g~^{-,t) depends continuously (in the 
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C"'^([0, D], L^(M))-topology) on the Gromov-Hausdorff distance. As the set 
M(A, Djig) is compact, the parameters C,N,I,K and 5 used in the above 
construction are uniform on M(A, D,io) and depend only on a. 

c. Construction of the boundary distance map. As 

Proposition m makes it possible to evaluate an approximate volume, vo1''(Mq). 
We will use this observation to construct a finite approximation to the set 
of the boundary distance functions associated with {M,g). Namely, for any 
X G M, let G L'^{dM) be given by 

rx{z) = dg{x,z), for any z G dM. 

The boundary distance functions determine the boundary distance map R : 

R(x) = r^. 

It turns out that the map R or, more precisely, its image R{M,g) is suffi- 
cient to reconstruct {M,g) (see [22] for the procedure of the reconstruction 
of (M, g) from R{M)). We will use an approximation R* to R{M, g) to con- 
struct the desired metric space (y, d) in Theorem ^. Note that the metric 
d is different to the induced metric as a (metric) subspace of L°°{dM) (cf. 
Step 1 in subsection d. 

To find R* consider subdomains C M, 

M; = {xeM: dg{x, Ti) G ((A - 2)r?, (A + 2)r?), / = 1, . . . , L} . (8) 

For any (3 the subdomain can be obtained as a finite number of unions, 
intersections and compliments of the subdomains Mq, of form (^. For any 

vol(J^'^) = vol(M) - vol(O), 

where QP = M \ and 

vol(17 n 0') = vol(17) vol(Sl') - vol(l^ U 0')- 

Moreover, for any a, /?, 

vo1(Mq, U M/j) = vol(M^), where 7; = max(a;, 

Therefore, it is possible to evaluate an approximate volume voP(M*), if we 
know {/Uj, V'j}- Clearly, |vol"(M^) - vol(M^)| < a when {^j, ^j} is 5-close 
to {Aj, ipj\dM} with sufficiently small 5. 
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When voP(M^) > cr we associate with this f3 a function rp G L°°{dM), 

rpiz) = Pir] for z gTi. 

Thus, if vol(M^) > 2a so that vol"(M^) > a, there is an approximate 
function rp such that for any x £ M^, 

\\rp -VxW < 2r]. 

In the construction of the approximate volume of a set ^ such that 
vol(M^) < 2a it may happen that vol"(Af^) < a. Thus, according to 
the procedure, R* do not contain rp corresponding to this multi-index /?. 
Consider next a point x G for such /?. By condition (^) there is an 
upper bound N{r]) for the number of all multi-indexes /?. Moreover, due 
to condition iv) of definition ^, vo\{BM{x,r])) > crj"^. Therefore, if u < 
C' if^ N {ri)~^ , there is a multiindex ^,^^13 such that 

i. vol(M*) > 2a, so that vor(M*) > a, 

ii. d{x,M*) < r]. 
Thus 

1^7 -rx\\ < 4??. 

Summarizing the previous considerations we obtain the following result. 

Lemma 1 For any (M, g) £ M. and any r/ > there is 5 > such that if 
{yUj, tpj} is 6-close to the boundary spectral data {Xj, ^jIom} then the set 
R* constructed from tf^j] j = 1, . . . , n{6~^)} is rj-close to R{M,g), i.e. 

dH{R*,R{M,g))<r,. 

Here dn is the Hausdorff distance between subsets in the metric space L°°{dM). 
Moreover, 6 can be chosen uniformly on any M(A,D,io)- 



d. Construction of an e-net The set R* has a metric structure inherited 
from L°°{dM), 



doo{ri,r2) = \\ri - r2\ 



L°°{dM)- 



If {M,g) is strongly geodesic (e.g. [IC], fS^]), i.e. all geodesic inside M are 
minimal, then 



(9) 
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for any x,y G M. Thus, the metric space {R*,doo) provides a 2?7-net for 

{M,g), 

dGH{iR*,d^), {M,g))<2rj. 

However, condition (^) is not vahd for general manifolds. Therefore, doo- 
metric is inappropriate to construct an approximation to (M, g) in the 
Gromov-Hausdorff topology. In this section we will describe how to equip 
the finite metric space R* with another metric d* so that {Y,d) = {R*,d*) 
becomes close to (M, g) in the Gromov-Hausdorff metric. In our description 
we will assume that the whole R{M) rather then R* is known. We will 
show how to find an approximation to dg{x,y) from R{M). This procedure 
is valid also in the case when we know only R* . The corresponding result is 
given at the end of this section. 

1. Assume first that xi,X2 G M*"* can be connected by a unique shortest 
geodesic, 7(3:1, 0:2) ■ Assume in addition that this geodesic can be continued 
beyond X2 to a boundary point z G dM so that 7(xi, z) is a shortest geodesic 
between xi and z. Then, 

dg{xi,X2) = dg{xi,z) - dg{x2,z) = \\rxi - r^all; 

can be found from R{M,g). 

It remains to identify those ri,r2 S R{M) which correspond to xi,X2 
with this property. As i*g on dM is known, it is possible to extend (M, g) 
to a larger Riemannian manifold {M,g) which has a continuous, piecewise 

metric g. We can then extend each function r^. G C{dM) to the function 
Tx G C(M \ M), rx{y) = dj^{x,y). Take a point z G dM and consider 
minimal geodesies j{xi,z), i = 1,2. Assume that 7(xj,z) \ {z} C M*"*, z 
is not a cut point along any of 7(xi, z) and 7(xj, z) are transversal to dM. 
Then these geodesies can be extended as minimal ^odesics beyond z upto 
some points yi G M \ M. Denote by ji{z,yi) C M \ M these extensions. 
Then the geodesic j{xi,X2) connecting xi and X2 extends to a minimal 
geodesic to z G dM if and only if 71 (z, yi) = 72(2, ^2)- 

2. Let now xi,X2 G M*"* lie on the sides 7(2/3,2/1) and 7(2/3,2/2) of a good 
geodesic triangle A (2/1, 2/2, 2/3)- A geodesic triangle A is called good if all 
three geodesies 7(2/2, 2/i) can be continued to dM as the shortest geodesies 
(see Fig. 2). Then, by part 1 we can find dg{yi,yj) and also dg{xi,yi). Let 
A(yi,y2,2/3) be a triangle in with de{yi,yj) = dg{yi,yj), where 4 stands 
for the Euclidean distance. Let Xj, i = 1, 2 be the points on the sides [y^, , yi] 
of A with de{xi,yi) = dg{xi,yi). Then, by the Alexandrov lemma, 

\dg{xi,X2) - 4(^1, ^2)1 < CCT^, 
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Figure 2: Example of a good triangle used to construct d{xi,X2)- In this 
case xi = Hi and X2 lies on the geodesic 7(^21^3)- 



where a = max dg{yi,yj) and c is uniform on M(A, D,io)- 

Therefore, if for some given xi,X2, we can find a proper geodesic triangle 
^(yi) 2/2, 2/3) of a small size, we can find dg{xi,X2) with an error of order cj^. 

3. The existence of a good geodesic triangle for any sufficiently close xi,X2 G 
M*"* is based on the following result. 

Lemma 2 Let A,D,io > and {M,g) e M{A,D,io). Let x £ M. Then 
there is a constant p > such that for any {M,g) £ M(A,D,io) cLnd x G 
{M,g) with dg{x,dM) > Cp there is a vector u £ T^M \u\ = 1 with the 
following property: 

Let {y,v) G Bp{x,u), \v\ = 1. Then the geodesic jy{tv), t > 0, can be 
extended as a shortest geodesic to the boundary. 

Here Bp{x,u) is the ball of radius p with center at (x,u) in the Sasakian 
metric on TM. Combining steps 1-3 we can find an approximation to 
dg{xi,X2) when xi,X2 are sufficiently close to each other and not too close 
to dM. (For the definition of the Sasakian metric, see e.g. [pS|].) 

4. Let 

dg{xi, dM) = min ri{z) < r] < io, ri = r^,- 

zddM 

There are unique Zi G dM with dg{xi,Zi) = ri{zi). Denote 

d = d{xi,X2) = [dlM{zi,Z2) + |ri(zi) - r2(z2)P]^^^ , 
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where dgMi', •) is the distance along dM. Then, 

\dg{xi,X2) - d{xi,X2)\ < C{r]d + d^). 

Therefore, d is a good approximation to d when Xi are close to dM. 

5. Steps 1-4 define approximate distance d{xi,X2) for sufficiently close 
xi, X2- Using a standard procedure we can extend d{xi,X2) onto arbitrary 
Xi G M. 

When we know the metric space {R*,doo) rather then R{M,g) we can 
carry out the above constructions approximately. This gives rise to the 
following result. 

Lemma 3 Let Y C L°°{dM) be a set described in Lemma |^. Lt is possible 
to equip Y with a metric structure d{ri,r2) so that 

dGH{{Y,d),{M,g))<e{ri), 

where e{rf) when — > 0. The function e{r]) can be found uniformly for 
M{A,D,io). 



Combining steps a.-d. we obtain a procedure to construct a finite e-net 
for {M,g),e = e{6). 

VII. Concluding remarks. 

Remark 8. The above considerations lack quantitative estimate for the 
function e(5). Combining our construction with the Carleman estimates 
(see, e.g. [^], ||3^, we will obtain quantitative estimates for e = £{6). 
This will be described in a forthcoming paper. 

Remark 9. All considerations remain valid for the Dirichlet Laplacian. 
We need just to take into account that 7^ in M*"*. 

Remark 10. It is possible to use the convergence of the heat kernels 
restricted to dM instead of the boundary spectral convergence. Besides, it 
is possible to modify the procedure to find an approximation to the heat 
kernel associated with — rather then R*. This puts our results within 
the framework of the spectral convergence of Riemannian manifolds (see e.g. 

[|l^). This will be described in the forthcoming paper. 

Remark 11. Requirements of C^-regularity of Rm and S are of tech- 
nical nature. We intend to improve reconstruction procedure to be valid for 
weaker regularity assumptions for Rm and S. 

Part of our results also briefiy announced in [20|. 
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